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Abstract 

The universal temperature dependence of metabolic rates has been used to predict how ocean biology will respond to 
ocean warming. Determining the temperature sensitivity of phytoplankton metabolism and growth is of special importance 
because this group of organisms is responsible for nearly half of global primary production, sustains most marine food 
webs, and contributes to regulate the exchange of CO2 between the ocean and the atmosphere. Phytoplankton growth 
rates increase with temperature under optimal growth conditions in the laboratory, but it is unclear whether the same 
degree of temperature dependence exists in nature, where resources are often limiting. Here we use concurrent 
measurements of phytoplankton biomass and carbon fixation rates in polar, temperate and tropical regions to determine 
the role of temperature and resource supply in controlling the large-scale variability of />i situ metabolic rates. We identify a 
biogeographic pattern in phytoplankton metabolic rates, which increase from the oligotrophic subtropical gyres to 
temperate regions and then coastal waters. Variability in phytoplankton growth is driven by changes in resource supply and 
appears to be independent of seawater temperature. The lack of temperature sensitivity of realized phytoplankton growth 
is consistent with the limited applicability of Arrhenius enzymatic kinetics when substrate concentrations are low. Our 
results suggest that, due to widespread resource limitation in the ocean, the direct effect of sea surface warming upon 
phytoplankton growth and productivity may be smaller than anticipated. 
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Introduction 

Temperature governs the metabolism of all organisms. Within 
favourable thermal ranges, there exists a positive, often exponen- 
tial, increase in metabolic rate as temperature rises [1,2]. On-going 
and predicted ocean warming [3] prompts the question of how 
phytoplankton photosynthetic activity and growth may respond to 
increasing sea surface temperatures [4—6]. This response will be 
relevant for the functioning of the Earth system as a whole, 
because phytoplankton sustain most marine food webs and, being 
responsible for nearly half of the global primary production [7] , 
contribute to control the exchange of CO2 and other radiatively 
active gases between the ocean and the atmosphere [8,9] . 

Laboratory work with cultures growing under optimal condi- 
tions has shown that the maximum growth rates of phytoplankton, 
and thus their maximum rate of biomass-specific production, 
increase exponentially with temperature with a Qjo of approxi- 
mately 2 within the tolerable temperature range [10-13]. 
However, it is not clear whether phytoplankton assemblages in 
nature show the same degree of temperature dependence in their 
realized rates of metabolic activity and growth. An analysis of vitro 
oxygen evolution data across the global ocean suggests that 
phytoplankton photosynthesis does increase with temperature and. 



furthermore, that the degree of temperature dependence is similar 
to that predicted by the metabolic theory of ecology [14]. 
Similarly, an eco-evolutionary model predicts that phytoplankton 
living in tropical regions sustain faster growth rates than their 
counterparts living in temperate and polar regions [6]. These 
studies, however, did not consider the role of resource supply, 
which is included in applications of the metabolic theory of 
ecology [1,15] and limits phytoplankton production and growth in 
most of the open ocean [16]. Nutrient limitation has been shown 
to reduce the temperature sensitivity of carbon fixation by 
phytoplankton [17-19]. If resource limitation attenuates the 
temperature dependence of metabolic rates, different regions 
where phytoplankton experience different degrees of resource 
limitation may show contrasting productivity responses to 
warming, and models that use laboratory-based values of Qjo 
are likely to overestimate the direct effects of temperature upon 
algal growth. 

Understanding the variability and controlling factors of 
phytoplankton growth rates is key to predict the response of 
ocean biology to external forcings such as climate variability [20] . 
However, in spite of decades of observation [10,20-23], a general 
picture of phytoplankton growth variation over broad spatial scales 
in the ocean remains elusive. Here, we use concurrent data of 
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phytoplankton carbon biomass and photosynthetic carbon fixa- 
tion, determined in coastal and open-ocean waters of polar, 
temperate and tropical biomes, to investigate the large-scale 

variability of phytoplankton metabolic rates, represented by the 
biomass-specific rate of carbon fixation [10,24]. Our analysis 
allows us to characterize biogeographic patterns in phytoplankton 
growth and to assess the relative role of temperature and resource 
availability in the control of marine primary productivity. 

Materials and Methods 

Data acquisition 

We compiled data on phytoplankton carbon biomass concen- 
tration and primary production rate that were determined 
concurrendy in surface assemblages of coastal and open-ocean 

regions. Coastal observations (number of sampling visits, n = 26) 
were conducted throughout the year at a central station 
(depth = 40 m) in the Ria de Vigo (NW Iberian Peninsula), a 
highly productive embayment subject to frequent upweUing 
events, particularly in spring and summer [25,26]. Open-ocean 
stations (re = 38) were visited in April-May 1996 and September- 
October 1996 along the Atlantic Meridional Transect, which 
crosses temperate, subtropical and tropical regions in the north 
and south Adantic Ocean [27]. The data from the AMT stations 
were grouped according to their latitude into north temperate (35- 
49°N), south temperate (35-48°S), north oligotrophic (20-3 TN), 
south oligotrophic (10-34°S) and equatorial and Mauritanian 
upweUing (5°S-2()°N). Additional data (re = 8) were obtained from 
the Soiree (Southern Ocean Iron Release Experiment) study, 
conducted in open-ocean, polar waters of the Australasian-Pacific 
Southern Ocean [28-30]. To minimise the effects of vertical 
variability in irradiance upon phytoplankton photosynthetic 
activity, the present analysis uses only data from surface samples 
(0-5 m in coastal waters, 0-20 m in open-ocean waters). 

Phytoplankton biomass and production 

Phytoplankton carbon biomass was estimated from measure- 
ments on fixed samples of cell abundance and biovolume, 
determined with flow cytometry for the picophytoplankton and 
with an inverted microscope for the nano- and micro-phytoplank- 
ton, following protocols that are detailed elsewhere [26,27,30]. 
Biovolume data were transformed into carbon biomass by 
applying appropriate conversion factors. Primary production rate 
was determined with the '*C -uptake technique using on-deck 
incubations which simulated in situ temperature and irradiance 
conditions [26,27,29]. Phytoplankton biomass turnover rates 
(d~'), equivalent to intrinsic growth rates [10,24], were calculated 
by dividing the daily carbon fixation rate by phytoplankton carbon 
biomass. 

Hydrography data 

Vertical profiles of temperature and salinity were obtained with 
a CTD probe attached to a rosette sampling system. Density was 
calculated from temperature and saKnity by using the standard 
UNESCO equation. The concentration of dissolved nitrate in the 
euphotic layer was determined by segmented-flow analysis using 
an automatic analyser and standard colorimetric protocols 
[26,27,31]. The detection limit of this analysis (0.05 flvcidi L~') 
was sufficient to measure nitrate concentration at the base of the 
euphotic layer, which was always >0.2 /flnol L The depth of 



the euphotic layer was determined from vertical profiles of 
irradiance obtained with a PAR Li-Cor sensor. 

Index of resource supply 

We computed a resource supply index {RSI) that takes into 
account the concentration of nitrate at the base of the euphotic 
zone (JV05[io/„p^-|), the density difference between the surface and 
the base of the euphotic zone (Act,), the depth of the euphotic 
zon(;, dcfinc-d as the 1 % PAR level (1 %PARz}, and the depth of the 
upper mixed layer [JJMLz), defined as the first depth at which ct, is 
0.125 units higher than the surface value [32]: 



The first term in Eq. 1 reflects the fact that nutrient transport to 
the upper layer increases with increasing nutrient concentration 
below, but also that this transport becomes progressively more 
limited as vertical stratification intensifies. The second term in Eq. 

1 ser\'es to differentiate conditions in which phytoplankton are 
confined to a shaUow upper mixed layer, relative to the depth of 
the euphotic zone, and thus experience relatively high average 
irradiances, from conditions in which intense mixing results in 
mixed layers whose depth equals or exceeds the penetration of 
irradiance into tlu; water column, pott'ntially resulting in light 
limitation of phytoplankton growth. Our index of resource supply 
is not intended to provide an accurate estimate of local-scale 
resource availability, but is used here to depict broad differences 
between regions with widely contrasting hydrographic regimes. 
RSIwsls not calculated for the stations sampled in the Soiree study, 
where, in addition to macronutrients and irradiance, iron was a 
limiting factor for phytoplankton production and growth [28] . 

Statistical analyses 

Pearson's r was calculated to assess the existence of linear 
dependence between variables. We applied reduced major-axis 
regression to determine the linear relationship bctw(X'n logig- 
transformed variables. Ninety-five percent confidence intervals for 
the regression parameters were calculated by bootstrapping over 
cases (2,000 repetitions). The non-parametric Mann-Whitney U 
test was used to compare means between groups of stations. 

Results 

Phytoplankton biomass and production 

Across all studied regions, which showed widely contrasting 
hydrographic properties (Table 1), carbon fixation rates covaried 
closely with phytoplankton standing stocks, represented either by 
carbon biomass (Fig. la) or by chlorophyll a concentration (chl a) 
(Fig. lb). The range of variability in phytoplankton production 
rates was >3 orders of magnitude, compared with approximately 

2 orders of magnitude for phytoplankton carbon biomass and chl 
a. The key feature of these relationships is that in both cases the 
slope of the fitted line is significantly higher than 1 , taking a value 
of 1.46 (95% confidence interval: 1.35, 1.58) in the production- 
carbon biomass regression and 1.29 (1.19, 1.39) in the production- 
chl a regression. This result implies that, as phytoplankton 
standing stocks increase, primary production rates increase- at a 
faster rate. From oligotrophic regions to coastal waters, mean 
phytoplankton chl a and biomass increased by a factor of 20-30, 
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Figure 1. Phytoplankton chlorophyll a, biomass and produc- 
tion. Scatterplots showing the log-log relationship between (a) primary 
production and phytoplankton carbon biomass, (b) primary production 
and chlorophyll a concentration and (c) phytoplankton carbon biomass 
and chlorophyll a concentration. Different symbols indicate the 
sampling region as defined in the IVlethods section. 'Soiree In' and 
'Soiree Out' refer to stations located inside and outside, respectively, 
the iron-fertilized patch during the Soiree study. Data were fitted to a 
linear model using reduced major axis regression. 
doi:10.1371/journal.pone.0099312.g001 



whereas primary production increased by a factor of >100 
(Table 2). 

The relationship between chl a and phytoplankton carbon 
(Fig. Ic) was remarkably robust, considering the uncertainties 
involved in estimating biomass from biovolume and the fact that 
the samples were analysed by several independent laboratories. 
The slope of the log-log regression between phytoplankton carbon 
and chl a (0.89; 95% CI: 0.83, 0.94) was significantly lower than 1, 
indicating that the C:Chl a ratio tends to decrease as phytoplank- 
ton standing stocks increase from open-ocean oligotrophic to 
open-ocean temperate and then coastal waters (Table 2). This 
pattern likely reflects the fact that in temperate and coastal waters 
phytoplankton are acclimated to lower mean irradiances, which 
results in increased chl a cellular content. 

Temperature and phytoplankton growth 

We found no overall relationship between phytoplankton 
growth rate, represented by the biomass turnover rate (P*"), and 
temperature (Pearson's ? = —0.045, p>0.3, n=12) (Fig. 2). The 
highest P*"' values (>1 d ') were measured in coastal waters, 
despite the fact that their temperature was 10-15°C lower than 
that of subtropical and tropical oligotrophic regions (Table 1). 
Some of the coldest and the warmest waters analysed had similarly 
low P*^' values (<0.4 d~ '). During the Soiree study, P*^' inside the 
iron-fertilised patch took higher values than in surrounding waters 
with the same temperature. The stations within the Atlantic 
upwelling region tended to have higher P*^ values than oligotro- 
phic stations with similar temperatures. 

Resource supply and phytoplankton growth 

P*"' increased with chl a (Fig. 3a), suggesting that an enhanced 
supply of resources, necessary to sustain larger standing stocks, 
leads to faster phytoplankton growth. Groups of samples with 
mean chl a above 2 mg m had mean P*^ values higher than 
1.2 d whereas samples with mean chl a below 0.2 mg m ^ 
showed P*^ values below 0.4 d~'. These differences were 
statistically significant (Mann- Whitney's U test, /)<0.05). The role 
of resource supply in coiitroUiiig phytoplankton biomass turnover 
was confirmed by the positive correlation between the resource 
supply index [RSI) and P*^ (Pearson's r = 0.45, /)<0.001, n = 59; 
Fig. 3b). RSI was calculated as a function of nitrate concentration 
at the base of the euphotic layer, the density difference between 
sea surface and the base of the euphotic layer, and the depths of 
the upper mixed layer and the euphotic zone (see Material and 
Methods). RSI varied widely, reflecting large across-system 
differences in nutrient concentration and water column stratifica- 
tion (Table 1). The oKgotrophic regions showed low i?>S7 values (< 
10 mmolN kg~') (Fig. 3c), mainly due to their low nitrate 
concentration at the base of the euphotic layer (Table 1). Higher 
nitrate concentration in the temperate regions resulted in RSI 
values above 20 mmolN kg~ . The highest RSI values were 
determined in the Ria de Vigo, where high nitrate concentrations 
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at the base of the euphotic layer co-occur with upper mixed layers 
~ ^ tliat are shallow relative to the depth of the euphotic layer 

■g_^ (Table 1). Despite having similar nitrate concentration at the base 

of the euphotic zone, the temperate regions had lower RSI values 
■B !| than the coastal waters (Fig. 3c), because the former were 

° s characterized by deeper mixed layers relative to the depth of the 

2 I euphotic zone. Finally, RSIwas lower in the upwelling region than 

,c o in the coastal waters of Ria de Vigo, even though deep nitrate 

concentration was higher in the former, because of the strong 
g I thermal stratification, associated with the upwelling, that resulted 

•t g in a relatively large density gradient (Table 1). 

2 I ^ There was an overall positive relationship between P"^' and RSI 

e ^ across all studied regions (Fig. 3c). The open-ocean temperate and 
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g S determined in the temperate and upwelling open-ocean regions 

and in the coastal waters (Mann-Whitney's U test, p<0.05). The 
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tions. Moreover, experimental studies with natural phytoplankton 
assemblages have shown that under nutrient limitation the 
temperature dependence of photosynthetic rate is reduced [17] 
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Table 2. Phytoplankton properties at the studied regions. 



Location 


Chi a 


Phyto C 


C:Chl a 


P 


Temperate N {n = 8) 


0.44±0.56 


30±22 


99 ±42 


13±15 


Temperate S {n-4) 


0.71 ±0.29 


34±10 


50±10 


15±6 


Temperate S {n = 4) 


0.71 ±0.29 


34±10 


50±10 


15±6 


Oligotrophic N (n = 5) 


0.10±0.03 


n±3 


117±15 


2.5 ±0.8 


Oligotrophic S (n = 13) 


0.13±0.07 


9±2 


83±31 


2.1 ±0.8 


Upwelling {n = 8) 


0.20±0.07 


16±11 


77±36 


10±11 


Coastal {n = 26) 


3.5±3.2 


202 ±187 


58±28 


248±271 


Soiree In (n = 6) 


1.19±0.16 


89±33 


75±25 


25±3 


Soiree Out (n = 2) 


0.20±0.03 


12±3 


65 ±27 


2.0 ±0.6 



Mean (± standard deviation) surface values of chlorophyll a concentration (Chi a, mg m"^), phytoplankton carbon biomass (Phyto C, mgC m"^), phytoplankton carbon 
to chlorophyll a ratio (C:Chl a), and rate of primary production (P, mgC m"^ d"^) in each location, n is the number of stations visited at each location. 
doi:10.1371/journal.pone.0099312.t002 



irradiance and increasing nutrient availability [45-47] . Therefore, 
normalising metabolic rates by chlorophyll a results in biased 
patterns of temperature dependence whenever there is covariation 
between temperature, nutrient availability and irradiance. 

Our observations depict a consistent, broad-scale pattern in 
phytoplankton growth rates, which increase from oligotrophic 
gyres to temperate, open-ocean regions and then to coastal, 
productive waters. This resource-driven pattern is associated with 
changes in community structure [26,27], as oligotrophic waters are 
dominated by small cells (picocyanobacteria and picoeukaryotes) 
with relatively low biomass turnover rates, whereas the most 
productive waters are dominated by fast-growing diatoms. Warm 
temperatures and a tight coupling between zooplankton-mediated 
nutrient regeneration and phytoplankton uptake could, in 
principle, result in an absence of nutrient limitation and near- 
maximal phytoplankton growth rates in the nutrient-impoverished 
subtropical gyres [20,21,23]. This possibility, however, contrasts 
with the experimental evidence showing physiological stress of the 
dominant picophotoautotrophs in the oligotrophic Atlantic, which 
is alleviated upon addition of nitrogen [16,48,49]. Our results 
support the existence of nutrient limitation of both standing stocks 
and growth rates in the oligotrophic gyres [50], i.e. a situation in 
which both Liebig's and Blackman's types of limitation apply 
[16,51]. This limitation appears to be less acute in temperate 
open-ocean regions, which are subject to seasonal nutrient 
injections through vertical mixing, and in the upwelling open- 
ocean region, where upward water motion enhances the nutrient 
supply to the euphotic zone. In the coastal, productive waters of 
Ria de Vigo, high resource supply results in growth rates that are, 
on average, similar to the maximum growth rates expected given 
the in situ temperatures. These resource-driven biogeographic 
patterns imply a predominantiy bottom-up control of phytoplank- 
ton growth in the ocean. 

A previous study [52] has shown that phytoplankton growth 
rates in coastal, nutrient-rich waters show a degree of temperature 
dependence comparable to that observed in laboratory cultures, 
which confirms that temperature sets an upper limit to the 
maximum growth rates that phytoplankton can achieve under 
resource-saturated conditions. In addition, nutrient stoichiometry 
and pigment content have been shown to be temperature- 
dependent, reflecting intracellular changes in the abundance of 
macromolecules involved in light harvesting, photochemistry and 

doi:10.1371/journal.pone.0099312.g002 



The lack of relationship between temperature and phytoplank- 
ton growth rate in our data may seem at odds with the results of 
some experimental studies in open-ocean waters which have 
shown that phytoplankton metabolism responds to temperature 
changes [42,43] . However, these manipulation experiments reflect 
short-term, transient responses in phytoplankton physiology and 
therefore do not capture the large acclimation potential of 
phytoplankton [44]. The present results are also in contrast to 
the findings of Regaudie-de-Gioux and Duarte [14], who 
concluded that there is a strong temperature dependence of 
phytoplankton photosynthesis in the sea. The discrepancy is likely 
to result from the fact that in their analysis these authors pooled 
together data obtained throughout the euphotic zone, thus 
introducing irradiance as a covarying factor. Samples near the 
base of the euphotic zone, where temperatures are colder than at 
the surface, will necessarily show lower rates of photosynthesis, not 
as a result of the lower temperature but as a result of lower 
irradiance. In addition, photosynthetic rates in the study by 
Regaudie-de-Gioux and Duarte were normalised by chlorophyll a 
instead of carbon biomass. However, it is well-established that the 
ceUular content of chlorophyll a increases with decreasing 
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Figure 2. Temperature and phytoplankton growth. Carbon- 
specific production rate (P"") plotted against temperature for all stations. 
Symbols as in Fiq. 2. 
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Figure 3. Resources and phytoplankton growth. Carbon-specific 
production rate (P*") as a function of (a) chlorophyll a concentration (chl 
a) and (b, c) resource supply index {RSI). P*" data in (a) were grouped 
into the following chl a ranges: 0-0.2, 0.2-0.5, 0.5-2, 2-5 and >5 mg 
m"^. Bars indicate ±1 standard deviation. Data groups labelled with 
different letters in (a) and (c) have significantly different P"" values 
(Mann-Whitney U test, p<0.05). The units of RSI derive from the division 
of nitrate concentration by seawater density. Note that RSI is not a 
nutrient concentration, but reflects both nutrient and light availability 
(see Material and Methods for details). 
doi:10.1371/journal.pone.0099312.g003 



biosynthesis [12,53,54]. Furthermore, inter-species differences in 
thermal tolerance ranges and optimal temperature for growth [55] 
contribute to explain changes in species distribution, community 
structure and associated ecosystem fianctioning over multiple 
temporal and spatial scales [56-59]. Nevertheless, as far as the 
large-scale variability in phytoplankton production and growth 
rates is concerned, temperature seems to be unimportant whereas 
resource supply plays a crucial role. This result may have 
important consequences for the modelling of ocean ecosystem 
responses to climate change. Biogeochemical models typically 
estimate phytoplankton growth rates as a multiplicative function of 
nutrient concentration, irradiance, and temperature-dependent 
maximum growth rate [5,60,61]. This formulation assumes that 
temperature effects on growth are independent of resource supply 
status, so that for a given condition of resource availability, an 
increase in temperature will necessarily result in faster phyto- 
plankton growth. In contrast, our observations indicate that 
resource limitation drastically reduces the temperature depen- 
dence of phytoplankton metabolism. Similarly, the growth [62] 
and extracellular enzymatic activities [63] of natural marine 
bacterial assemblages show only a very limited response to 
temperature when resources are in low supply. Hence, the 
sensitivity of microbial metabolic rates to sea surface warming, 
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Figure 4. Temperature and phytoplankton growth in samples 
with different resource availability. Carbon-specific production 
rate (P"") is plotted against temperature for oligotrophic stations with 
low resource supply index (fiS/<10 mmolN kg"') and temperate, 
upwelling and coastal stations with intermediate RSI (10-20 mmolN 
kg"'). Both linear regression fits were non-significant (p>0.6). 
doi:10.1371/journal.pone.0099312.g004 
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particularly in open-ocean, nutrient-limited regions, may be 
smaller than anticipated [4,14,64,65]. 
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